Honeycomb antiferromagnet with a triply degenerate dimer ground state 



On 
O 
O 

?3 



a 

o 
o 



> 

00 

in 
<N 
o 

o 
a^ 
o 



X 



Rakesh Kumar,* Dushyant Kumar, and Brijesh Kumar 
School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India 

(Dated: June 1, 2009) 

We present an antiferromagnetic quantum spin- 1/2 model on honeycomb lattice. It has two parts, 
one of which is the usual nearest-neighbor Heisenberg model. The other part is a certain multiple 
spin interaction term, introduced by us, which is exactly solvable for the ground state. Without 
the Heisenberg part, the model has an exact threefold degenerate dimer ground state. This exact 
ground state is also noted to exist for the general spin-S* case. For the spin-1/2 case, we further 
carry out the triplon analysis in the ground state, to study the competition between the Heisenberg 
and the multiple spin interactions. This approximate calculation exhibits a continuous quantum 
phase transition from the dimer order to Neel order. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.50.Ee, 75.40.Mg 



I. INTRODUCTION 

The low dimensional quantum spin systems are a sub- 
ject of great current interest. Much of this research di- 
rectly attends to the real (quasi) one and two dimensional 
spin systems studied in the laboratories^'^'^'^. There is 
also a formal side to it which is concerned with investigat- 
ing, at various theoretical levels, the effects of low spatial 
dimensionality, quantum spin fluctuations and frustra- 
tion on the nature of the ground state of a model spin 
system. Since the antiferromagnetism is sensitive to all 
of these, the ground state of a quantum antiferromagnet 
(AF) can choose from a variety of possibilities (known, or 
not yet known) For example, the ground state of the 
spin-1/2 nearest neighbor (nn) Heisenberg antiferromag- 
net is a critical spin-liquid (with power law decay of the 
spin-spin correlations and zero local magnetic moment) 
in one dimension (Id), but it has Neel order on two and 
higher dimensional bipartite lattices. Furthermore, the 
competing interactions can induce changes in the nature 
of ground state, say, from being a Neel ordered state to 
becoming spontaneously dimerized^ (or something else). 
The Majumdar- Ghosh model presents an exactly solvable 
case of a spontaneously dimerized doubly degenerate sin- 
glet ground state in Id^. Similar spin models have also 
been constructed in two dimensions (2d)^'-^^'-^-^'-^^'-^^. 

Of the spin systems in 2d, the honeycomb lattice comes 
across as a special case to study. Its site-coordination 
is three which lies between Id and the square lattice. 
Hence, the quantum fluctuations are expected to be 
stronger on the honeycomb than on the square lattice. 
Moreover, the honeycomb is not a br avals lattice. It has 
two spins per unit cell. Thus, we expect a natural case for 
spontaneous dimer order, without breaking translational 
symmetry, in a honeycomb antiferromagnet. For the 
spin-1/2 nn Heisenberg AF on honeycomb lattice, various 
calculations^^' ^^'^^ indeed show larger quantum fluctua- 
tions than on the square lattice, but the ground state still 
exhibits Neel order (although weaker than square lattice). 
Several other studies have shown that under various frus- 
trated conditions the ground state on honeycomb lattice 
can get disordered^^'^^'^^'^^'^^. Motivated by these ob- 



servations, one of us (BK) constructed a quantum spin 
model with multiple spin interactions on honeycomb lat- 
tice, which has an exact triply degenerate dimer ground 
state (see Fig. 2). Here, we present this model, and inves- 
tigate using triplon mean-fleld theory the transition from 
the dimer to Neel order in the ground state, in the pres- 
ence of the nn Heisenberg AF interaction. This model 
could be of further interest in investigating deconflned 
quantum criticality in 2d antiferromagnets^'^^'^^. 

The paper is organized as follows. In Sec. II, we discuss 
the model and its exact ground state. In Sec. Ill, we do 
the triplon mean-fleld theory. In Sec. IV, we discuss the 
results. Finally, we conclude with a summary. 



II. MODEL 

In this paper, we study the following quantum spin-1/2 
model on honeycomb lattice (pictorially shown in Fig. 1). 
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O 
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Pl2^34^56 



■ ^23^45^61] (1) 



Here, Sf^- = (S^ + Sj)^. The flrst term above is the near- 
est neighbor Heisenberg model with interaction J. The 




FIG. 1: (Color online) The red and green bullets represent the 
spins on two sublattices of the honeycomb lattice. The con- 
necting lines denote the nearest neighbor antiferromagnetic 
interaction, J. The shaded hexagon denotes the multiple spin 
interaction, present on every hexagonal plaquette. 
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second term is what we have introduced to reahze the 
dimer ground state (as in Fig. 2). This multiple spin 
interaction is generated by the product of the pairwise 
total spins of three nn pairs on each hexagon [in two 
different ways, (1, 2)(3, 4)(5, 6) or (2, 3)(4, 5)(6, 1)]. It in- 
volves all the six spins of an hexagon, and the summa- 
tion is taken over all hexagonal plaquettes of the hon- 
eycomb lattice with periodic boundary condition. Af- 
ter expanding Sf^- as | + 28^ • Sj (for a pair of spin- 
1/2), and regrouping different terms, Eq. (1) becomes: 
H = ^KL + ^ ^(4) ^ ^(6)^ ^Yieie L is the total 

number of lattice sites, and H^'^\ H^^^ and H^^^ denote 
the quadratic, quart ic, and sextic spin interactions, re- 
spectively. Note that H respects the lattice translation 
and point-group symmetries, and is also SU (2) invariant. 
We take J^K > 0, and J -\- K = 1 sets the unit of energy. 
Thus, J =1-K, where K e [0, 1]. 



/ 



/ 
/ 
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A. Exact ground state 



Consider the model for = 1. It can be written as: 

(2) 



= ^yZ i^Kl + hK2] 

^ o 



S?9SLSL and hK2 



SiaSlsSii. In this 



where Hri - 012^34^56 
case, the nn Heisenberg exchange is absent. We only 
have the multiple spin interactions. Clearly, hxi and 
hK2 have positive eigenvalues, with zero as the minimum. 
Hence, the ground state energy of Hk is bounded below 
by zero. The operator hxi gives zero when at least one 
of three concerned spin pairs, that is (1,2), (3,4) or (5,6), 
forms a singlet. Similarly for hK2' Therefore, on a single 
hexagon, a zero energy eigenstate of hKi-\-hK2 can be ob- 
tained by simultaneously forming singlets on the opposite 
edges. It leaves the remaining two spins remain as 'free'. 
For example, one such state is [1, 2] Ims) (g) [4, 5] Ime), 
where [ij] = (| tiij) - | ii^j))/V^ is the singlet formed 
by i and j spins, and rrik =T or |. Moreover, there are 
three ways of choosing such dimer forming spin pairs. 
After knowing these single plaquette dimer states, it is 
straightforward to show that the three dimer ordered 
configurations shown in Fig. 2 form the exact zero en- 
ergy ground state of Hk on the full lattice. We have 
cross-checked it using numerical diagonalization on a fi- 
nite spin cluster in the following subsection. These dimer 
states have been known to arise in the ground state of 
the quantum dimer model on honeycomb lattice, but ours 
is probably the first example of a SU (2) spin model on 
honeycomb with this dimer ground state^^. 

Let the ground state configurations be denoted as 
102), and 1^3). These dimer states do not break the trans- 
lational symmetry of the lattice, and are obviously SU (2) 
invariant. The point group rotational symmetry is bro- 
ken, however. The wave function of the dimer state, 
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FIG. 2: (Color online) The dimer states |(/)i), 102), and |03) 
form the ground state of the Hamiltonian, Hk- The thick 
blue lines denote the dimer singlets. 



can be explicitly written as: 



101 



n [^'^']' 



(3) 



where V is the set of singlet forming dimers in the state 
|0i). The other two states, 102) and |03), are related to 
|0i) via the threefold rotation as: 



1^3) 



C3 |0l) 

Cl \4>i) 



(4a) 
(4b) 



where C3 is the clockwise 27r/3 rotation operator. 

To this end, we would like to mention that this exact 
ground state of Hk is also valid for the general spin- 
S system. For the spin-^* case, a dimer would denote 
a singlet state formed by a pair of spin-^*. Everything 
else (that is, the dimer pattern, the degeneracy, and the 
ground state energy) is the same. Since the maximum 
total spin of a pair of spin-S* is 25*, we can rescale the 
coupling ^ to 2(25+1^(25+1) ' ^his just makes the energy 
contribution of the multiple spin interaction comparable 
(in powers of S) to that of the Heisenberg part. 
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B. Finite size numerical diagonalization 

We have done exact numerical diagonalization of H 
on a 12-site honeycomb cluster with periodic boundary 
conditions. This is just to numerically verify the exact 
ground state on a small cluster. Our present expertise 
does not allow us to do exact numerical diagonalization 
on larger spin clusters. We only use total magnetization 
and spin-inversion symmetries in the coding. The ex- 
act diagonalization results clearly show that the ground 
state for the exactly solvable case {K = 1) is indeed triply 
degenerate with zero energy (see Fig. 3). The next eigen- 
state has a finite energy gap to the ground state. Away 
from the exact case, the ground state energy decreases 
smoothly without level crossing. Although the degener- 
acy seems to lift as soon as K is different from 1, but we 
expect the degeneracy to survive, in a finite range of K 
values, for large enough systems. 

0.2, ^ , ^ , ^ , ^ , ^ , 



LU 




K 



FIG. 3: (Color online) Energy eigenvalues per site from the 
exact diagonalization of if on a 12-site honeycomb cluster. 

In order to ascertain the nature oi K = 1 ground state, 
we compute the spin-spin correlation, (S^ • Sj), and the 
four-point (dimer-dimer) correlation, 

D(i,j,k,i) = ((Si ■ Si)(Sfe ■ SO) - (Si ■ S,-){Sfc ■ S,). (5) 

The latter helps in identifying the dimer order. Accord- 
ing to this definition, the Dij^k,i is positive when the 
two dimers are correlated, and negative if the dimers 
are uncorrelated. We compute these first in the numer- 
ically generated ground state, and then compare them 
with those calculated in the exactly known dimer ground 
state. On the infinite lattice, the nn spin-spin correlation 
is equal to —1/4 (further neighbor spin correlations are 
identically zero in the states of Fig. 2), and the dimer- 
dimer correlation is 1/8 when two dimers are perfect sin- 
glets. 

For = 1, the numerical ground state wavefunctions 
would be some orthogonal linear combinations of 
1^2), and 1^3). The choices of the linear combination 
are not unique, however. As we have not implemented 
translation and point group symmetries in our compu- 
tational scheme, this ambiguity in the degenerate out- 
put states of our (less sophisticated) program remains. 



TABLE I: Spin-spin correlations 



Dimer 


Numerical diagonalization 


Exact 


(2,1) 


-0.244768 


-0.244768 


(2,3) 


-0.245889 


-0.245889 


(2,4) 


-0.008221 


-0.008221 


(2,5) 


-0.248879 


-0.248879 


(2,6) 


-U.UUzz4z 


-U.UUzz4z 


(2,7) 


0.000373 


0.000373 


(2,8) 


-0.000373 


-0.000373 


(2,9) 


0.002242 


0.002242 


(2,10) 


-0.002242 


-0.002242 


(2,11) 


0.000373 


0.000373 


(2,12) 


-0.000373 


-0.000373 


TABLE II: Dimer-dimer correlations. 


Dimers 


Numerical diagonalization 


Exact 


(2,5)(3,4) 


-0.062599 


-0.062225 


(2,5)(4,1) 


-0.061476 


-0.061414 


(2,5)(6,7) 


-0.061476 


-0.061414 


(2,5)(7,8) 


-0.062599 


-0.062225 


(2,5)(9,10) 


-0.062599 


-0.062225 


(2,5)(10,11) 


-0.061476 


-0.061414 


(2,5)(11,12) 


-0.062599 


-0.062225 


(2,5)(12,9) 


-0.061476 


-0.061414 


(2,5)(4,7) 


0.122756 


0.124064 


(2,5)(6,9) 


0.124718 


0.124718 


(2,5)(8,11) 


0.124718 


0.124718 


(2,5)(10,1) 


0.124718 


0.124718 


(2,5)(12,3) 


0.124718 


0.124718 



Therefore, we use the zero temperature thermal density 
operator to correctly compute the ground state proper- 
ties. For an operator O, its thermal average is given by 
(O) = Tr(pO), where p = Z'^e'^^ is the thermal den- 
sity operator {Z = Tr e~^^). In the zero temperature 
limit, the density operator reduces to 

1 

'^=]^El*'^)(*''l (6) 

where Ng is the degeneracy of the ground state, and l^^y) 
are the ortho-normalized ground state eigenvectors. In 
the present calculation, Ng = 3. 

Since the numerical eigenstates are orthonormal, we 
use them directly to compute the correlations in the 
ground state, as prescribed above. This data is shown 
in the second column of the Tables I and II. The exact 
wavefunctions, etc. are not orthogonal. Therefore, 
we first orthogonalize them using Gram-Schmidt proce- 
dure (on the same cluster as used for numerical diago- 
nalization; see Fig. 4), then apply the density operator 
averaging to compute the correlations. These are given 
in the third column of the two tables. Clearly, the spin 
correlations are nearest neighbor type, and the dimer cor- 
relation matches with dimer order in the exact ground 
state (compare Fig. 2 with Fig. 4). Numbers from the 
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exact and the numerical calculations match perfectly. 




FIG. 4: (Color online) The dimer-dimer correlations are cal- 
culated by taking the dimer (2,5) as reference. Thickness of 
the dimers is proportional to their dimer correlation values. 
The blue dimer denotes the positive correlation (i.e., a sin- 
glet), whereas purple ones represent negative correlation. 



III. TRIPLON MEAN-FIELD THEORY 

While at = 1, the ground state of H has an exact 
dimer order, but it is known to be a Neel ordered AF 
state when K = {). It would be interesting, therefore, 
to make some investigation of the transition from the 
dimer to Neel ordered ground state, as K is varied. Here, 
we present an approximate study of this quantum phase 
transition by doing triplon analysis with respect to the 
dimer phase. A triplon is a triplet excitation residing 
on a dimer, and dispersing according to the interactions 
present in the system. While a non-zero gap in the triplon 
dispersion corresponds to dimer phase, the gaplessness 



implies AF order in the ground state^'^^. 

The triplon analysis is conveniently carried out in the 
bond operator representation in which the singlet and 
three triplet states of a pair of spin- 1/2 (a bond) are 
described in terms of the corresponding bosons (called 
bond operators)^^. The bosonic creation operators, s"^ 
and t]^ {a = x, ^, z), respectively create singlet or triplet 
states on a bond, subjected to the physical constraint 
^'''5 + tl^ta = 1 (repeated Greek indices summed over). 
The two spins on a dimer are represented as 

Sia = ^ (^sHa + 4s - ieaf3jt^ptj^ (7a) 
S2a = ^ (^-sha - t^S - iCafd-ftpt-f^ (7b) 

where subscripts 1 and 2 denote, say, left and right spins, 
and Cafd-f is the totally antisymmetric tensor. In the sim- 
plest triplon analysis, the singlet background is treated 
as a mean field {s^ = 5), and the triplon dispersion is cal- 
culated by ignoring the triplon-triplon interaction (please 
find a better description of the triplon analysis in Ref. 8; 
we follow the same strategy as therein.). 

Although the exact dimer ground state is triply degen- 
erate, we can only take one of these as the reference state 
to do the triplon mean-field calculation. We rewrite H in 
terms of the bond operators, taking as the reference 
dimer background. The Heisenberg exchange on a dimer 
can now be written as: 

Si(R).S2(R) = -^s-2 + itti^, (8) 

where R is the position vector of the dimer. For the spins 
coming from different dimers, we have 



S^(R)-S,(R + 5) = (-l)^+'^ j[ti(R)t^{IL + 6_) + ta{R.)ta{R + 5.) + h.c.] p = l,2). 



(9) 



By using Eqs. (8) and (9), and applying the constraint on 
bond operators globally (with /i as Lagrange multiplier), 
we get the following mean-field Hamiltonian for H. 



Hmf = i?0 + i5]{(A-S-2a) 



ka'^-ko 



.Jkc]} (10) 



and 



cos I -/c.^ I COS I —^k, 



y I ' 



(12) 



(13) 



Here, the triplon operators have been Fourier trans- 
formed from the dimer lattice to the corresponding re- 
ciprocal lattice, with k vectors lying in its first Brillouin diagonal mean-field Hamiltonian can be written as 
zone. Moreover, 



where A is the effective chemical potential. Equation (10) 
is diagonalized by the Bogoliubov transformation. The 
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l + -K-{J+-K]s'--X + \s 



, (11) 



Hmf = £^0 + ^ -Ek ( 7L7ka 



(14) 



5 



where 7kQ;'s are Bogoliubov bo sons, and the triplon dis- 
persion, = Y^A (A — > 0. The ground state 
energy per site is given by 



eo 



_3_ 
2L 



(15) 



where eo = Eq/L. By minimizing the ground state en- 
ergy with respect to A and s^, we get the self-consistent 
equations, whose solution gives the mean-field results. 



A. Dimer phase 

When the minimum of the triplon dispersion is 
nonzero, the dimer phase is stable against triplet excita- 
tions. Therefore, the gapped triplon phase corresponds 
to having dimer ground state. The self-consistent equa- 
tions is this case are: 



k 

9 „ 3A ^ r?k 



-K 



k 



(16a) 
(16b) 



where r^k = — ^Ks'^ cos (^^kx) cos (^^ky^ . These 

equations are obtained by minimizing the ground state 
energy, i.e., dcg/dX = and deg/ds^ = 0. 

The weight of having singlet state on a dimer is mea- 
sured by . If all the dimers form perfect singlets (like in 
the exact case), then = 1. Otherwise, we get 5^ < 1, 
due to triplon fluctuations in the ground state. 



B. Neel phase 

As K is gradually decreased away from if = 1, at some 
point we find that the triplon gap vanishes (see Fig. 5). 
The triplon dispersion now touches zero at k = Q, where 
Q = (0, 0). It means the triplon occupancy at wavevector 
Q becomes singular, which implies the Bose condensation 
of triplons at Q. Thus, we need to introduce a third 
quantity (in addition A and 5^), the triplon condensate 
density, nc, which is notionally given by 



2 3 ^A-^^^Q 



-c = ^(t^.tQ.).-^^-^J. (17) 
Now, the revised set of self-consistent equations are: 



1 

2^ 



3A 
L 



L.^_ Eu 



5 o 

k#Q 



(18a) 
(18b) 

(18c) 



Physically, the non-zero ric corresponds to AF order in 
the ground state, which in the present case comes out to 
be the Neel order. The staggered magnetic moment in 
the Neel phase is given by, Ms = s^/n^. 



IV. RESULTS AND DISCUSSION 

Now, we present the results obtained by the triplon 
mean- field calculation. The self-consistent Eqs. (16) of 
the gapped phase are solved for A and 5^, for different 
values of K. Interestingly, for = 1, it gives = 1, 
same as the exact answer (see Fig. 8). At the exact point, 
the triplon dispersion, E'k (plotted in Fig. 5), is flat with 
an energy gap of 1.125. This value of gap at = 1 
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FIG. 5: (Color online) Triplon dispersion Ek for different val- 
ues of K. The r, A, and B denote the wavevectors (0,0), 
(27r/3,0), and (27r/3, 27r/3V3) in the Brillouin zone. 

is in agreement with a direct estimate of 9/8, which is 
calculated as: A(x=i) = {^\Hk\^) / {^\^) ^ where 




{k,l}. 



(19) 



Here, V is a set of all dimers in V except (/c, /), and {/c, /} 
denotes a triplet state on dimer (/c, /). The set V is same 
as in [see Eq. (3)]. 

For K < 1, £^k acquires a finite width, with minimum 
at the F point [that is, Q]. Eventually, it touches zero 
at Q, and remains so, below K* = 0.256. For K < K* , 
we solve the Eqs. (18). As shown in Fig. 6, now the 
staggered magnetization acquires a non-zero value while 
the gap remains zero. 

The triplon mean-field theory thus predicts a contin- 
uous transition from the dimer to Neel ordered phase in 
the ground state of H for spin- 1/2. The phase diagram 
is just a line presented in Fig. 7. As pointed out earlier, 
the K = 1 model has the same exact ground state for 
higher spins also. Therefore, in some future studies, it 
would be interesting to extend this quantum phase dia- 
gram to include a spin- axis, with S = 1/2, 1,3/2,... to 
S ^ oo (the classical limit). The classical case can be 
discussed right away. Let the spins be classical vectors. 
The Hk will now have an infinitude of spin configurations 
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FIG. 6: (Color online) The spin gap, A = ^q, and the stag- 
gered moment Mg. Inset: the triplon condensate density, Uc. 




FIG. 9: (Color online) The ground state energy per site from 
the mean-field calculation. It compares reasonably with the 
exact diagonalization calculation (see Fig. 3). 
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FIG. 7: Mean-field quantum phase diagram of H for spin-1/2. 



in the ground state (not related via global spin rotation), 
because the spins on each dimer (of Fig. 2), separately, 
must cancel. Hence, Hk itself is a frustrated model. But 
it does not compete against the classical nn Heisenberg 
interaction for winning the ground state as the infinite 
set of pairwise spin-cancelled configurations also includes 
the Neel states. Therefore, i^* = 1 in the classical limit. 
This discussion reveals an important feature of H that 
is, the multiple spin and Heisenberg interactions don't 
compete against (or frustrate) each other. Instead, the 
quantum mechanics acts better when K is sufficiently 
large. Hence, the Neel to dimer transition in the ground 
state of H is driven purely by quantum fluctuations. 




FIG. 8: (Color online) The singlet weight, s^, and the effective 
chemical potential A. 

Below, we present the mean-field critical behavior of 
the spin-gap and staggered magnetization. The critical 



exponent for both is 1/2, which is derived by analyzing 
the gap and Uc equations in the small neighborhood of 



A ^ 0.122(i^-i^*)^/2 
Ms ^ 0.44(i^* - i^)^/2 



(20a) 
(20b) 



Fig. 10 shows an enlarged plot of the mean-field data 
around i^*, together with the estimated critical behavior. 
The two compare well. 



0.03 



0.02 



0.01 - 



0.00 



--- Ms 

• Ms-0.440 (K*-K)1/2 

— A 

■ 0.122 (K-K*)1/2 



"ooo" 



0.01 

(K-K*) 



0.02 



FIG. 10: (Color online) The spin-gap and staggered mag- 
netization near the critical point. The estimated results of 
Eqs. (20) are also plotted. 



V. SUMMARY 

We have constructed and studied a quantum spin-1/2 
model on honeycomb lattice. In one limit of the interac- 
tion parameter {K =1), the model has an exact threefold 
degenerate dimer ground state. Away from the exact case 
(0 < < 1), we study the evolution of the ground state 
using triplon mean-field theory. The mean-field theory is 
exact at K = 1 , and it shows a continuous quantum phase 
transition from the dimer-ordered to Neel ordered ground 
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state at iC* = 0.256. Within this mean-field theory, the 
critical exponents for the spin-gap (in the dimer phase) 
and the staggered magnetization (in the Neel phase) are 
1/2. We have also done preliminary numerical calcula- 
tions on this model. We have done exact diagonalization 
on a 12-site honeycomb cluster. For = 1, it gives 
the triply degenerate ground state with the same spin- 
spin and dimer-dimer correlations as in the exact dimer 
ground state. We are now focusing on doing numerical 



work on larger spin clusters. 
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